How many bits does it take to track an open quantum system? 
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A D-dimensional Markovian open quantum system will undergo quantum jumps between pure 
states, if we can monitor the bath to which it is coupled with sufficient precision. In general these 
jumps, plus the between-jump evolution, create a trajectory which passes through infinitely many 
different pure states. Here we show that, for any ergodic master equation, one can expect to find 
an adaptive monitoring scheme on the bath that can confine the system state to jumping between 
only K states, for some K > (D — l) 2 + 1. For D = 2 we explicitly construct a 2-state ensemble for 
any ergodic master equation, showing that one bit is always sufficient to track a qubit. 

PACS numbers: 03.65.Yz, 03.65.Aa, 42.50.Lc, 42.50.Dv 



The first quantitative model of quantum dynamics was 
Einstein's model of stimulated and spontaneous jumps [lj 
between Bohr's stationary atomic states [H. In modern 
language, this is a model for an open quantum system 
weakly coupled to a heat bath, and, as Einstein showed, 
such jumps can lead to an equilibrium state that is a 
thermal mixture of energy eigenstates, with finite en- 
tropy. In Einstein's model, if one could track the individ- 
ual stochastic events of energy exchange between atom 
and bath, then one would know which energy eigenstate 
the system occupied at any time. Truncating to a finite 
number D of energy eigenstates, it follows that only a 
finite classical memory is required to keep track of the 
quantum system (that is, to know its exact pure state) 
in thermal equilibrium: a K -state memory with K = D. 

Einstein's theory is a special case of Markovian 
open quantum system dynamics for finite-dimensional 
systems, which most generally are describable by a 
Lindblad-form master equation (ME) 
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Cp = -i[H eS p - pHlft] + ^2 cipc], 



(1) 



where H e s = H — i J^i cJc;/2. Here H is Hermitian (it is 
the Hamiltonian) but the jump operators {£;} are com- 
pletely arbitrary. Einstein's theory is a special case be- 
cause in it each jump operator is proportional to \E) (E'\, 
for some -ff-eigenstates \E) and \E'), so that the state 
after any jump is a stationary state \E). For a general 
ME, it is always possible, in principle, to monitor the 
bath such that every jump is resolvable, so that the sys- 
tem can be known to be in some pure state \ip{t)) a t 
all times 0-0] ■ However, in general, after a jump at 
time Tj, the state cx q \ip(Tj)) will depend on the prc- 
jump state \ijj(Tj)), and will not be an an eigenstate of 
H. Even if it were an energy eigenstate, it would not in 
general remain stationary until the next jump, because 
its subsequent evolution would be generated by the ef- 
fective (non-Hermitian) Hamiltonian H e ff appearing in 
Eq. (I]). 

It is thus not at all obvious whether for a general finite- 



dimensional open quantum system it would be possible to 
keep track of its pure state, even in principle, with a finite 
classical memory. On the face of it, it would seem neces- 
sary to store the exact times of each jump — a sequence 
of real numbers {tj : j} each of which would require, in 
principle, an infinite memory to store. Alternately one 
could store the conditioned quantum state \ip{t)} itself, 
but this (a D-dimensional complex vector) would also 
require an infinite memory. This situation is of course 
completely different from a finite-state stochastic classi- 
cal system (which is what Einstein's model amounts to), 
where a finite-state classical memory of the same size is 
always sufficient. 

In this Letter we address this fundamental question 
about open quantum system dynamics. We show that for 
any ergodic Markovian dynamics [15J of a D-dimensional 
quantum system, one can expect to be able to track 
the state with a if -state classical apparatus for some 
K > (D — l) 2 + 1. This is possible only because there is 
entanglement between the system and bath, which means 
that different monitoring schemes on the bath give rise to 
different sorts of stochastic pure state trajectories ( "un- 
ravellings" for a given ME. We then prove that for 
D = 2 (a qubit), K = 2 is always sufficient; that is, 
there is always an unravelling for which the qubit jumps 
between only two possible states, \4>\) and \(j>2). Al- 
though this sounds similar to Einstein's dynamics, it is 
in fact quite different in general — the two states are 
non- orthogonal, (4>\\4>2) ^ 0, and the monitoring of the 
qubit 's environment must be adaptive, controlled by the 
classical bit that stores the state of the qubit. 

We begin by revisiting the preferred ensemble fact [||, 
to explain why it is not possible in general to unravel a 
ME such that the system jumps between the eigenstates 
of the equilibrium p (as in Einstein's model). Then we 
show the general result cited above for D-dimensional 
systems, and give an explicit construction of the adap- 
tive unravelling for the special case of cyclic jumps with 
only one Lindblad operator c. We then prove that one 
bit is always sufficient to track a qubit, and that in some 
cases it is actually possible to store the state of the open 
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qubit using less than one bit of memory on average, if 
we imagine an ensemble of N 3> 1 qubits, each indepen- 
dently measured. Surprisingly, considering K > 2 can 
actually help in this regard. We illustrate these phenom- 
ena using the resonance fluorescence ME 

The preferred ensemble fact. - Consider a Linbla- 
dian £ with unique steady state defined by Cp ss = 
and we assume that p ss is a mixed state. This mixed 
state can be decomposed in terms of pure states \<f>f.) 
via p ss = EfcLi Pk \4>k) (4>k\ with positive constants p k . 
Note that there are infinitely many such decompositions, 
as the states \cj>k) need not be orthogonal. However, only 
for some decompositions is it possible to devise a way 
to monitor the system's environment — which leaves the 
average evolution of the system unchanged from Eq. (fTJ) 
- such that the system will only ever be in one of the 
states \4>k)i an( i w iU spend a proportion of time in that 
state equal to p k in the long-time limit. Decompositions 
{pk, \<Pk)} that can be realized in this way are called phys- 
ically realizable (PR). The fact that some decompositions 
are not PR is known as the preferred ensemble fact 

As shown in Ref. |f|, an ensemble {p k , | is PR iff 
(if and only if) there exists rates Kj k > such that 

K 

Vfc, C\cf> k ) (<f> k \ = 5> ifc (l&> - l&> (M) ■ (2) 
fc=l 

For a general ME, most decompositions {pk, \<Pk)} of p ss 
arc not PR, including the K = D ensemble composed 
from the diagonal basis for p ss Q . 

The existence of PR ensembles. - - For finite K and 
D, searching for solutions of Eq. (|2|) reduces to solving 
polynomial equations. We can describe K pure states 
with K(2D — 1) real unknowns and K quadratic con- 
straints from normalization. Eq. ([2]) introduces K 2 — K 
unknown rates Kjk and imposes an additional K(D 2 — 1) 
cubic constraints (the minus one is because both sides 
are traceless by construction). Thus we have KD 2 poly- 
nomial constraints and K(2D + K — 2) unknowns. For 
K > (D — l) 2 + 1 we have an underdetermined system of 
equations. For linear underdetermined systems, given by 
equations {fj}, there are infinitely many solutions except 
for the set of measure zero for which there exist constants 
{otj} such that y] ■ ocjfj+1 = 0. Similarly, for polynomial 
systems {pj}, the Real Nullstellensatz certifies that 
there are no real solutions iff there exist some polynomi- 
als {aj} and {d k } such that yj • a,jPj + yj fc d 2 , + 1 = 0. 
Thus, for a general ME we expect to be able to find a 
A'-element PR ensemble for some K — 1 > (D — l) 2 . 

The freedom that experimentalists have (in princi- 
ple) to realize different PR ensembles comes from the 
ability to monitor the system's environment in differ- 
ent ways. This can be understood as follows. The ME 
(TT]) is invariant under the transformations {£;} — » {c^}, 
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— ft m c' m '). Here ft is an arbitrary com- 



plex vector and S is an arbitrary semi-unitary matrix — 
E™=i Sp m S m i = Sv.i. Unravelling this ME with {c' m } 
as the jump operators and H eS = H' — i E m =i Kn&mfe 
as the effective non-Hcrmitian Hamiltonian clearly gives 
different stochastic evolution, while leaving the average 
evolution unchanged. To obtain the most general pure- 
state unravelling of the ME, we require ft and S to de- 
pend upon the previous record of jumps. That is, we 
require an adaptive monitoring 0, [H| . Of course when 
we use this to achieve jumping between a finite number 
of states, the classical iC-state memory that stores which 
state the system is currently in carries all the information 
necessary for determining ft and S. That is, the adaptive 
unravelling is specified by K different values for ft and S. 
The physical meaning of these parameters is most eas- 
ily explained in a quantum optics context: S describes 
a linear interferometer [161 ] taking the field outputs from 
the system as inputs, while ft describes adding (weak) 
local oscillators to the output fields from the interfer- 
ometer prior to detection by photon counting. Recently, 
adaptive control of a weak local oscillator has been used 
for optimally distinguishing coherent states Q , and of a 
strong local oscillator for improved phase estimation 0] . 

Backing out the measurement scheme. - Although 
for every PR ensemble there must exist a monitoring 
scheme exists by definition, it may not be easy to find. 
We now present an explicit method for determining this 
scheme for the special case of cyclic jumps with a single 
Lindblad operator c. That is, we assume that the sys- 
tem in the state \<j>k) always jumps to the state \4>k+i) 
(strictly, |0(fc+i) modif ))• I n this case there are only K 
jump rates, so the number of real unknowns is only 2KD. 
With a single Linblad operator and cyclic jumps, both 
sides of Eq. @ have rank two by construction, so for 
D > 2 Eq. @ is less constraining than in the general 
case. Nevertheless, the system will be overconstrained 
for D > 2. Thus in general we do not expect there to 
exist cyclic jump solutions for D > 2. Later we exhibit 
eleven different cyclic jump solutions for a qubit (D = 2), 
and the method here is applicable to each of them. 

In the case of a single Lindblad operator, the only free- 
dom in the unravelling is in choosing ft k , the local oscil- 
lator amplitude when the system is known to be in state 
\<j>k) (that is, when the if -state classical memory is in 
state k) . This gives the jump operator c + ft k and the ef- 



fective Hamiltonian H„ 



H eS + i, 
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l k \ 2 /2. Thus 



the system will undergo cyclic jumps iff 



where 
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{H eS + ift k *c)\<t> k ) cx \<f> k ) , (c + ft k )\0 k ) cx \</> k+1 ), (3) 

Now by assumption from Eq. ([2]), C \<j>k) {4>k \ oc 
\<pk+i) (4>k+i \ — \4>k) (4>k\- From this we can show that 

c\<j) k ) = a k \<j> k ) +b k \<j> k+1 ) , (4) 
H c ft \<f>k) = Cfc \4> k ) + ia* k b k \<j>k+\) ■ (5) 
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For some coefficients a,k, b k and c k that are easily found 
given \4>k) and \4>k+i). Comparing this to Eq. (|3]), we see 
that choosing f3 k = —a k gives cyclic jumps as required. 

For a qubit, one bit is all it takes. - - We now prove 
that a 2-state PR ensemble always exists for a qubit. We 
use the Bloch representation, so that Eq. (J]) becomes 



'= Af + b, 



(6) 



where A is a 3 x 3-matrix and b is a 3-vector. As al- 
ways, we assume that there exists a unique steady state 
f ss = —A^^^b, which is the case iff the real part of each 
eigenvalue of A is negative. We can track this sys- 
tem with a if -state memory iff there exists an ensemble 
{pfc, rfc} and rates Kjk > such that 



Vfc, f k ■ r k =1, 



A' 



Vj, Afj +b=^n jk (f k -fj). 



(7) 
(8) 



fe=i 



Thus the problem reduces to finding a real solution to 
a system of quadratic equations with real coefficients. 
This type of problem is surprisingly hard even for a small 
number of unknowns, and is known to be an NP-complcte 
problem in general |l0|. 

Luckily, the simplest case of K = 2 has an analytical 
solution. Here the qubit is assumed to jump between 
two states. r\ and 7" 2 . Then Eqs. (J8j) reduce to a single 
equation, A(f 1 - f 2 ) = [n\2 + K 2 i)(f 2 — ^1), which is 
simply an eigenvalue equation. Thus f\ — F2 = v is an 
eigenvector of A with eigenvalue A = — (K12 + k 2 i). From 
Eq. ([7]) and the fact that pi = K21/ («2i + ^12) = 1 — P2, 
it is simple to show that f\ and r 2 have the form 



ri =f ss + vy/1- ||r ss ||Vpi/(l-pi) (9) 
r 2 =r ss - vy/l- ||r ss ||V(l-pi)/pi (10) 



with v = v/\\v\\ is the normalized eigenvector, and 



Pi 



i/l - ||f ss || 2 + (r ss ,v)< 



(11) 



Because the Bloch vectors must be real, only real eigen- 
vectors v of A can contribute to the solution. By as- 
sumption, A has three nonzero eigenvalues and, by the 
fundamental theorem of algebra, at least one eigenvalue 
(and consequently one eigenvector) is real. Therefore, a 
qubit always has a preferred ensemble comprising just 
two states. 

Entropy. - As nooted above, it may be possible to 
store the state of a qubit in less than one bit, in an av- 
erage sense. We can quantify this by using the Shan- 
non entropy. Under continuous monitoring, in the long- 
time limit, the system will occupy states \<f>k) with prob- 
abilities pfe. The Shannon entropy for this ensemble is 



h ({pfc}) = - Efc Pfe log 2 Pk- This is lower bounded by 
the von Neumann entropy for the steady-state mixture: 



h({pk})>S( Ps , 



-Tr[p ss log 2 p s; 



(12) 



with equality iff {pfc, \4>k) (</>fc|} is the diagonal ensemble. 

Note that if the eigenvector v of matrix A used to con- 
struct the ensemble in Eqs. ([ MTU]) is orthogonal to the 
steady state r ss , then the probability of occupying states 
r*i and ?i is 1/2. In this case the Shannon entropy is 1, 
meaning that one bit is sufficient to track the state of the 
system. If, on the other hand, v is not orthogonal to f ss , 
then the Shannon entropy h for the ensemble will be less 
than one and one could store the state on the qubit in less 
than one bit on average. That is, one could keep track of 
the state of a collection of N identically monitored qubits 
using only Nh bits, in the limit of large N. 

Resonance fluorescence. - In order to illustrate our 
ideas, we consider the example of resonance fluoresence 
of a two- level atom (a qubit) with basis states |0) and |1). 
The atom is coupled to the continuum of electromagnetic 
radiation and so decays to |0) at rate 7. At the same time, 
it is driven by a classical field with Rabi frequency f2. The 
qubit evolution in the interaction frame is given by a ME 
of the form of Eq. {T} with H = fi(|0) (1| + |1) (OQ/2 and 
one jump operator c = ^/7|0) (1| [3j. For this case, A 

and b in the Bloch vector equation, Eq. (|6|) , are 
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and b 



(13) 



The steady state f ss = (0, 27ft, ~7 2 ) T / (j 2 + 2Q 2 ) , is 
a mixed state for f2 ^ 0. The eigenvectors of A arc 
Vi = (1, 00) T and v ± = (0,7 ± ^ 7 2 - 16^ 2 ,4fi) T . 
Eigenvector v\ is real and orthogonal to f ss . Thus it 
will always yield a solution with Shannon entropy h = 1. 
This solution was originally discovered in Ref . [7| • Eigen- 
vectors v± are real only for |f2| < 7/4 and yield solutions 
with h < 1. In fact, for e = tt 2 /j 2 <C 1, the entropy of the 
solution due to u_ differs from S(p ss ) only at 0(e 3 ). Thus 
for small driving, this ensemble has an entropy very close 
to the bound (|T2|) . This is seen in Fig. []] where dashed 
lines show h for the three different 2-state solutions. 

For e > 0.0625 (that is, |fi| > 7/4) there are no low- 
entropy 2-state solutions. Surprisingly, by increasing K , 
the number of states, from 2 to 3, we regain a relatively 
low entropy solution for some range of e > 0.0625. Recall 
that for 2-state jumping, we use one real eigenvector of A 
to construct the PR ensemble, as the Bloch vectors must 
be co-linear with r ss . For 3-state jumping, the Bloch vec- 
tors must be coplanar (and not colinear) with f ss , so we 
require two eigenvectors. When A has complex eigenvec- 
tors, they come in conjugate pairs, and we can, for some 
values of fi, construct a PR ensembles with 3 real Bloch 
vectors using these conjugate pairs. 
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FIG. 1: The average number of bits required to keep track of 
the pure state of a qubit described by Eq. (|13|) . as a function 
of dimensionless driving power e = Q 2 /j 2 . The solid line is 
the von Neumann entropy for p ss , a lower bound on the mem- 
ory required. The other nine curves are for eleven different 
adaptive unravellings. Dashed lines are for 2-state jumping. 
Dotted lines are for cyclic 3-state jumping. 

We find all possible 3-state cycles by numerical search 
for solutions to Eqs. ([7][5]) using symbolic-numerical algo- 
rithms based on computing Groebner bases (l2j . There 
are eight solutions, coming in pairs, as shown (as dot- 
ted lines) in Fig. Q] For e < 0.0795, there are two PR 
ensembles generated from complex eigenvectors v± of A. 
These give the highest and lowest of the 3-state jumping 
entropy curves. In the region e < 0.0610, there are an ad- 
ditional four solutions constructed from v\ and V- . Only 
two new curves appear in Fig. [1] because they come in 
degenerate pairs. Finally, for e < 0.0335, there are two 
more solutions constructed from v±. As |0| decreases, 
the entropy for four of the solutions approaches 1.206, 
whereas the entropy for the other four approaches 0. 

The three 2-state PR ensembles and two of the 3-state 
PR ensembles are shown in Fig. [2] Figure [2ja) shows the 
main features of low-entropy solutions: the states are far 
apart, and the qubit spends most of the time in one state 
that is nearly aligned with the steady state. Figure [5Jb) 
captures the nature of high-entropy solutions: the states 
cluster around steady state. In both of these cases the 
ensemble lies in the x = plane. The ensemble with 
h = 1, Fig. I^c) does not, but is still symmetric under 
reflection in this plane. This symmetry is respected for all 
ensembles found except for the pairs of 3-state jumping 
solutions that are degenerate with respect to entropy. 
Figure [2]Jd) shows one such pair: the two ensembles are 
mirror images of each-other in the x — plane. 

In summary, we considered an arbitrary ergodic 
Markovian open quantum system subject to continuous 
monitoring that resolves every jump and allows the sys- 
tem to stay in a pure state. Under a generic monitoring 
scheme the system state will explore a manifold of pure 
states, so tracking it would require infinite memory. Here 
we showed that this situation is not an intrinsic prop- 
erty of open quantum systems, but is just a consequence 




FIG. 2: Solid arrows show Bloch vectors for 2-state jumping 
(a)-(c), and 3-state jumping (d). Colors [online] match those 
of Fig. [T] The volume of the sphere at the tip of each arrow 
represents the probability that the qubit occupies the corre- 
sponding pure state. The dashed 2UTOW IS 7"ss . For all plots, 
f2/7 = 0.2 (e = 0.04). 

of using the "wrong" monitoring scheme — a finite (in- 
state) classical memory is sufficient to track the state 
of the system, by adaptivcly changing the scheme used 
to monitor the environment, controlled by the state of 
the classical memory that stores the state of the quan- 
tum system. In general one would expected to need at 
least K — (D — l) 2 + 1 classical states to track a D- 
dimensional quantum system. The gap between K and D 
may be related to the recent result that there are stochas- 
tic processes that can be generated using quantum sys- 
tems of lower dimensionality than is possible using only 
classical systems The above quadratic difference, 

K — 1 = (D — l) 2 , is also reminiscent of other com- 
parisons between quantum and classical systems (l4j . so 
whether this K is always sufficient is an important open 
question. For D = 2, however, the answer is now known, 
and is: yes, one bit is always enough to track the state 
of a qubit . 
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